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Using 385 fb _1 of e + e~ collision data collected at center-of-mass energies around 10.6 GeV, we 
search for time-integrated CP violation in the Cabibbo-suppressed decays D°/D° — > 7r~7r + 7r° and 
D / D° — > K~ K + ty° with both model-independent and model-dependent methods. Measurements 
of the asymmetries in amplitudes of flavor states and CP eigenstates provide constraints on theories 
beyond the standard model, some of which predict CP violation in amplitudes at the 1% level or 
higher. We find no evidence of CP violation and hence no conflict with the standard model. 
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PACS numbers: 14.40.Lb, 13.25.Ft, 11.30.Er 



Charge-parity violation (CPV) [l], manifested as an 
asymmetry between the decay rates of a particle and 
its CP-conjugate antiparticle, requires at least two in- 
terfering complex quantum mechanical amplitudes with 
different phases. The strong phase of each amplitude 
respects CP symmetry while the weak phase changes 
sign under charge-conjugation. In the standard model 
(SM), direct CPV is due to relative weak phases that 
typically enter as a difference in phase between "tree 
level" and "penguin" 0] SM amplitudes. The pen- 
guin amplitudes in charm decays are, however, too 
small (0(0.1%) 0) to provide significant CPV. Ex- 
tensions of the SM introduce additional amplitudes of 
0(1%) 0, 0, with relative weak phases that can pro- 
duce CPV in charmed particle decays [f|. Current ex- 
perimental searches 0, i, & m, el m arc approaching 

this level of sensitivity. Observation of CPV with current 
experimental sensitivities would provide strong evidence 
of new physics. 

A recent theory paper Q argues that singly Cabibbo- 
suppressed (SCS) D (meaning either D° or D°) decays 
are uniquely sensitive to CPV in c udd, uss transitions 
and probe contributions from supersymmetric gluonic 
penguins. Such transitions do not affect the Cabibbo- 
favored (c — > sdu) or doubly Cabibbo-suppressed (c — > 
dsu) decays. Time-integrated CP asymmetries in D de- 
cays can have three components: direct CPV in decays 
to specific states, indirect CPV in D°-D° mixing, and 
indirect CPV in interference of decays with and without 
mixing. Indirect CPV is predicted to be universal for 
amplitudes with final CP cigenstates, but direct CPV 
can be non-universal depending on the specifics of the 
new physics. 

We search for time-integrated CPV in the three-body 



SCS decays D 



~tt + tt°,R R + tt°. These decays pro- 



ceed via CP eigenstates (e. 



P°7T , 



bir°) and also via 
flavor states (e.g., p ± 7r T , R* ± R T ), thus making it pos- 
sible to probe CPV in both types of amplitudes and in 
the interference between them. Measuring interference 
effects in a Dalitz plot (DP) probes asymmetries in both 
the magnitudes and phases of the amplitudes, not sim- 
ply in the overall decay rates. We adopt four approaches 
in our search for evidence of CPV, three of which are 
model-independent. First, we quantify differences be- 
tween the D° and D° DPs in two dimensions. Second, 
we look for differences in the angular moments of the 
D° and D° intensity distributions. Third, in a model- 
dependent approach, we look for CPV in the amplitudes 
describing intermediate states in the D° and D° decays. 
Finally, we look for a phase-space-integrated asymmetry. 
The first two methods are sensitive to differences in the 
shapes of the D° and D° DPs, allowing regions of phase 
space with CPV to be identified. The third method as- 



sociates any CPV observed using the first two methods 
with specific intermediate amplitudes. The last method 
is insensitive to differences in the DP shapes, so comple- 
ments the other methods. To minimize bias, we finalize 
the analysis procedure without looking at the data. 

We perform the present analysis using 385 fb _1 of e + e _ 
collision data collected at 10.58 GcV and 10.54 GeV 
center-of-mass (CM) energies with the BABAR detec- 
tor [l3[ at the PEP-II storage rings. The event se- 
lection criteria are those used in our measurement of 
the branching ratios of the decays D — > 7r~7r + 7r° and 
D — > K~K + tt° [3]. In particular, we study D mesons 
produced in D* + — > L>°tt+ and D*~ — > D°Tr~ de- 
cays that distinguish between D° and D . We require 
the D candidate CM momentum > 2.77 GeV/c and 
\m D *± - m D - 145.4 McV/c 2 | < 0.6 MeV/c 2 . Here, 
m refers to a reconstructed invariant mass. Around 
±1 standard deviation of the nominal D mass, we find 
82468 ± 321 tt-tt+tt and 11278 ± 110 R-R+tt signal 
events with purities of about 98%. We determine the sig- 
nal reconstruction efficiency as a function of the position 
in the DP using simulated D° and D° decays [14[ from 
e + e~ — ► cc events, subjected to the same selection pro- 
cedure that is applied to the data. 

A direct comparison of the efficiency-corrected and 
background-subtracted DPs for D° and D° events is the 
simplest way to look for CPV. Figure [1] shows the nor- 
malized residuals A in DP area elements, where 



(n D0 - R ■ n D o) I ^jol_ a + 



R 2 



(1) 



and n denotes the number of events in a DP clement and 
a its uncertainty. The factor R, equal to 0.983 ±0.006 for 
tt-tt+tt and 1.020 ± 0.016 for R-R+tt , is the ratio of 
the number of efficiency-corrected D° to D° events. This 
is introduced to allow for any asymmetry in the produc- 
tion cross section due to higher order QED corrections 
or in the branching fractions for D° and D° decay to the 
same final state. 

We calculate % 2 jv = (J2i=i A 2 )/^, where v is the 
number of DP elements: 1429 for 7r~7r + 7r° and 726 for 
K~ K + tt°. In an ensemble of simulated experiments with 
no CPV, we find the distribution of x 2 l v values to have 
a mean of 1.012 ± 0.001 (1.021 ± 0.002) and an r.m.s. of 
0.018 (0.036) for tt-tt+tt (R'R+tt ). The measured 
value in the data is 1.020 for Tr-7r + Tr° and 1.056 for 
K~K +, n a , so we obtain a one-sided Gaussian confidence 
level (CL) for consistency with no CPV of 32.8% for 
7t-7t + 7t° and 16.6% for K~K + tt . The same analysis 
procedure, when applied to simulated samples with ei- 
ther 1% fractional change in magnitude or 1° change 
in phase between the D° and D° amplitudes for decay 
to any of the main resonant states, gives a x 2 l v that is 
about 2(7 away from the no CPV hypothesis. Systematic 
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m 2 (7r-7i°) (GeV 2 /c 4 ) 



m 2 (K"7C°)"(GeV 2 /c 4 ) 



FIG. 1: (color). Normalized residuals in Dalitz plot elements, defined in Eq. [TJ for (a) D — > ty 7r + 7r° and (b) D — > K K + tt°. 



uncertainties are small (as will be clear from the model- 
dependent results of Tables IT1 ITT|) and have not been in- 
cluded in the CL calculation. 

The angular moments of the cosine of the helicity an- 
gle of the D decay products reflect the spin and mass 
structure of intermediate resonant and nonresonant am- 
plitudes [Hj]. We define the helicity angle Oh for decays 
of the type D — > r{AB) C as the angle between the mo- 
mentum of A in the AB rest frame and the direction op- 
posite to the D momentum in that same frame. The an- 
gular moments (l6l | of order I are defined as the efficiency- 
corrected invariant mass distributions of events weighted 
by spherical harmonics Y^Qh) = yl/27r Pi(cos9h)- 
Here Pi are the Legendre polynomials of order To 
study differences between the D° and D° amplitudes, we 
calculate the quantities X/ for I = — 7, where 



correlation coefficient between X,-, X, 



X 



(Pi-R-Pi) 
-7T + i?2 -^ 



(2) 



and Pi (Pi) are obtained from D° (D°) events. Higher 
moments are zero within errors in both data and simu- 
lation. For illustration, we show the Xi distributions for 
1 = 0-2, in Fig. H 

We then define x 2 jv of the angular moment distribu- 
tions of a two-body channel summed over all intervals in 
invariant mass as 

\ i=0 j=0 



Pij 



{XiXi) - (X) (X,) 



(xf) - <x 4 } 2 • ^(x]) - ( Xj y 



(4) 



We determine the in each mass interval by simulat- 
ing experiments with no CPV . We test the method on 
real data by randomly assigning events as D Q or D a , and 
then calculating \ 2 l v for the difference in their angular 
moments. We repeat this experiment 500 times and find 
the resulting \ 2 l v distribution to be consistent with no 
CPV, validating our calculation of p^. We then look at 
the D flavor in the data and calculate the X^/v values 
for the two-body channels with charge combinations +, — 
and +,0. Finally, we obtain a one-sided Gaussian CL 
for consistency with no CPV using the reference value 
and r.m.s. deviation from simulation. We find the CL 
for no CPV to be 28.2% for the tt+tt", 28.4% for the 
tt+tt , 63.1% for the K+K~, and 23.8% for the K+n 
sub-systems. Again, a 1% fractional change in magni- 
tude or 1° change in phase of any of the main resonant 
amplitudes gives a x 2 /^ that is about 2a away from the 
no CPV hypothesis. 

The Dalitz plot amplitude A can be parametrized as 
a sum of amplitudes A r (s+, s_) for all relevant interme- 
diate states r, each with a complex coefficient, i.e., A = 
Yl r ar e ^ r A r (s+, S-), where a r and <fi r are real. Here 
s+ and S- are the squared invariant masses of the pair 
of final state particles with charge combinations +, and 
— ,0. The fit fraction for each process r is defined as 



f r = J \a r A r \ ds+ds-/ f \A\ ds + ds. 

CP-symmetric background empirically 15 



. We model incoherent, 
17||. In the 



where v = 8k, k is the number of intervals, and is the absence of CPV, we expect the values of a r and 4> r (and 
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FIG. 2: (color online). Normalized residuals for the first three Legendre polynomial moments of the ty~ty + (row 1), 7r + Tr° (row 
2), K~ K+ (row 3), and K + it° (row 4) sub-systems. The confidence level for no CP violation (dashed line) is obtained from 
the first eight moments. The error bars represent ±1<t. 



hence f r ) to be identical for D° and D° decay. The re- 
sults obtained with this assumption are listed in Ref. [l~7| 
for D -> tt'tt+tt and in Ref. for D -> R-R+tt . 
To allow the possibility of CPy in the present analysis, 
we let a second process - not necessarily of SM origin 
- contribute to each of the amplitudes A r , thus permit- 
ting the a r , (j> r , f r for D° and D° to differ. We sum- 
marize the results of_the fit to the data in terms of the 
differences Aa r = a® — , A<fi r = </> D — 4>® , and 
A/r = f?° - fr° in Table^for tt-tt+tt and in Tabled 
for R-R+tt . The CP asymmetry in any amplitude, rel- 
ative to that of the whole decay, is no larger than a few 
percent. 

Systematic uncertainties in the quantities describ- 
ing CP asymmetries, reported in Tables HTHTl arise from 
experimental effects, and also from uncertainties in the 
models used to describe the data. We determine these 
separately, as described in Refs. ljj 17 1, and add them 
in quadrature. For all variations described below, we 
assign the maximum deviation from the central value 
as a systematic uncertainty, accounting for correlations 
among parameters. For resonance lineshapes and form- 
factors, we vary the parameters fl8j by ±lc Similarly, 
we vary the signal efficiency parameters for separately 
for D° and D° events by ±lcr, the ratios of particle- 
identification rates in data and simulation by Ha, and 
the background shapes by using simulation rather than 
data sidebands. We include uncertainties from D°-D° 
misidentification, estimated from simulation, in the ex- 
perimental systematic uncertainty. 

To this point, we have described the investigation of 
time-integrated CP asymmetry in neutral D meson de- 
cays using information from the DP distributions. Dif- 
ferences in the overall branching fractions for the D° 
and D° decays to 7r~7r + 7r°, R-R+tt would also indi- 
cate time-integrated CPV. This information is not cap- 
tured by the differential comparisons of the DP struc- 



tures already described, and is complementary to them. 
To correct for any production asymmetry in D-flavor 
assignment, we weight each event by the relative effi- 
ciency for flavor assignment, as described in Ref. 0]. 
Since there is an asymmetry Q between the number of 
events reconstructed at forward and backward polar an- 
gles (0?jo) °f the D candidate CM momentum, we extract 
the CP asymmetry value, a CP = ^^S__bS ; \ n intervals of 
| cos 9 C pa | . Here, TV denotes the number of signal events. 
Any forward-backward asymmetry is canceled by aver- 
aging over symmetric intervals in cos O'ffi , as shown in 
Eqs. 3-5 of Ref. 0]. In Fig. [3] we show the a CP for events 
in the D mass window used in the DP analysis. We 
perform % 2 minimization to obtain the central values: 
[-0.31 ± 0.41 (stat) ± 0.17 (syst)] % for tt-tt+tt and 
[1.00 ± 1.67 (stat) ± 0.25 (syst)] % for R-R+tt final 
states. The systematic uncertainties result from signal 
efficiency, particle-identification, background treatment, 
and D°—D° misidentification. As a consistency check, we 
repeat the analysis with a larger D mass window (±2.5<r) 
and find consistent results: [-0.28 ± 0.34 (stat) ± 0.19 
(syst)] % for tt-tt+tt and [0.62 ± 1.24 (stat) ± 0.28 
(syst)] % for R-R+tt . 

In summary, our model-independent and model- 
dependent analyses show no evidence of CPV in the SCS 
decays D — > Tr-Tr+Tr and D — > R-R+tt . The interme- 
diate amplitudes include well-defined flavor states (e.g., 
P ± tt+, R*+R+) and CP-odd eigenstates (e^a., p°tt°, 
<j)TT°). With the null results of Ref. @, S, ©, E3| for CP- 
even eigenstates D — > R+R~ and D — > tt+tt - , we con- 
clude that any CPV in the SCS charm decays occurs at 
a rate which is not larger than a few percent. These re- 
sults are in accord with the SM predictions, and provide 
constraints on some models beyond the SM 0]. 

We are grateful for the excellent luminosity and ma- 
chine conditions provided by our PEP-II colleagues, and 
for the substantial dedicated effort from the comput- 
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FIG. 3: (color online). Phase-space-integrated CP asymmetry as a function of the cosine of the polar angle of the reconstructed 
D candidate CM momentum for (a) D — > tt~ty + tt° and (b) D — > K~ K + ir° decays. The dashed lines represent the central 
values, and the shaded regions the ler intervals. 



TABLE I: Model-dependent CP asymmetry in the D — > 7r~7r + 7r° Dalitz plots. The first and second errors are statistical and 
systematic, respectively. For details on the Dalitz plot parametrization and the a r , <f> r , and f r values, see Ref. [l7| . As explained 
in text, Af r is closely related to Aa r and A<j!> r 



State 


fr (%) 


Aa r (%) 


A<j> r (°) 


Af T (%) 


P + (770) 


68 


-3.2±1.7±0.8 


-0.8±1.0±1.0 


-1.6±1.1±0.4 


p°(770) 


26 


2.1±0.9±0.5 


0.8±1.0±0.4 


1.6±1.4±0.6 


p-(770) 


35 


2.0±1.1±0.8 


-0.6±0.9±0.4 


0.7±1.1±0.5 


p+(1450) 


0.1 


2±11±8 


-30±25±9 


0.0±0.1±0.1 


p°(1450) 


0.3 


13±8±6 


-1±14±3 


0.1±0.2±0.1 


p~(1450) 


1.8 


-3±6±5 


8±7±3 


-0.2±0.3±0.1 


p+(1700) 


4 


19±27±9 


9±7±3 


0.4±1.0±0.4 


p°(1700) 


5 


-31±20±12 


-7±6±2 


-1.3±0.8±0.3 


p~(1700) 


3 


-3±14±11 


-3±8±3 


-0.5±0.6±0.3 


/o(980) 


0.2 


0.0±0.1±0.2 


-3±7±4 


0.0±0.1±0.1 


/o(1370) 


0.4 


-0.3±1.3±1.2 


7±14±5 


-0.2±0.1±0.1 


/o(1500) 


0.4 


0.4±1.1±0.7 


-1±12±1 


0.0±0.1±0.1 


/o(1710) 


0.3 


-3±3±2 


-25±13±11 


0.0±0.1±0.1 


/ 2 (1270) 


1.3 


8±4±5 


2±5±2 


0.1±0.1±0.1 


ct(400) 


0.8 


-0.3±0.7±2.0 


-4±7±3 


-0.1±0.1±0.1 


Nonres 


0.8 


12±7±8 


11±9±4 


0.2±0.3±0.2 
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